arXiv:1503.08256vl [cond-mat.stat-mech] 28 Mar 2015 


EPJ manuscript No. 

(will be inserted by the editor) 


Deconfined criticality for the two-dimensional quantum 
S = 1-spin model with the three-spin and biquadratic interactions 

Yoshihiro Nishiyama 

Department of Physics, Faculty of Science, Okayama University, Okayama 700-8530, Japan 

Received: date / Revised version: date 

Abstract. The criticality between the nematic and valence-bond-solid (VBS) phases was investigated for 
the two-dimensional quantum S — 1-spin model with the three-spin and biquadratic interactions by means 
of the numerical diagonalization method. It is expected that the criticality belongs to a novel universality 
class, the so-called deconfined criticality, accompanied with unconventional critical indices. In this paper, we 
incorporate the three-spin interaction, and adjust the (redundant) interaction parameter so as to optimize 
the finite-size behavior. Treating the finite-size cluster with N < 20 spins, we estimate the correlation- 
length critical exponent as — 0.88(3). 

PACS. 75.10.Jm Quantized spin models - 05.30.-d Quantum statistical mechanics - 75.40.Mg Numerical 

nally, the idea was developed in the context of the 

gauge-field-theoretical description for the two-dimensional 
strongly-correlated-electron system. Meanwhile, it turned 
out that the underlying physics is common to a variety 
of systems in terms of the emergent gauge field 

As a lattice realization of the deconfined criticality, 
the quantum S = 1/2-spin square-lattice antiferromag- 
net with the plaquette four-spin interaction, the so-called 
J-Q model PlIS] . has been investigated extensively; the 


simulation studies - 74.25.Ha Magnetic properties 


1 Introduction 


The phase transition between the Neel and valence-bond- 
solid (VBS) phases for the two-dimensional quantum spin 
system is attracting much attention recently [Tl[2l|3l|4j[5l 

It is expected that the phase transition, the so-called de¬ 
confined criticality, belongs to a novel universality class, 
accompanied with unconventional critical indices. Origi- 
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bipartite-lattice systems such as the square- [HIS] and honeycombTo be specific, we present the Hamiltonian for the two- 


lattice antiferromagnets do not conflict with the 

quantum Monte Carlo method, and large-scale-simulation 
results are available. However, it is still unclear whether 
the phase transition is critical or a 

weak first-order transition with a latent-heat release [3 
[Tll^l^ fTnilTT] ■ The controversy may be reconciled by a re¬ 
cent renormalization-group analysis Ea, which revealed 
an influence of a notorious marginal operator around the 
deconfined-critical fixed point; see Ref. [53] as well. Be¬ 
cause the character of the marginal operator depends on 
each lattice realization, it may be sensible to survey a va¬ 
riety of lattice realizations. 

The S = 1-spin model is a clue to the realization of the 
deconfined criticality UMsm. A key ingredient is that 
the S = 1-spin model admits the biquadratic interaction, 
which stabilizes the VBS phase as the spatial anisotropy 
varies m- The phase transition separating the VBS and 
nematic [53] phases is expected to belong to the deconfined 
criticality [5^ . We consider a non-bipartite-lattice version 
(Fig. [T]); the details and underlying ideas are explained 
afterward. In the preceding paper |26| . the correlation- 
length critical exponent was estimated as v = 0.92(10). In 
this paper, based on the preceding studies [541I551I56] . we 
incorporate a rather novel interaction term intrinsic to the 
S = 1-spin model, namely, the three-spin-interaction term 
[5i|3D] (in addition to the biquadratic one), and survey the 
extended parameter space so as to optimize the finite-size 
behavior. 


dimensional S = 1-spin model 

n = . s, + (s,. s,)2] -J'Y. (s,. s,)2 

{ij) <(b)) 

+J"{ ^ S,. S, 

{(ij)) 

+j3j2[{S.-S,){S,-S,) + h.c.]}. ( 1 ) 

[ijk] 

Here, the symbol Si denotes the quantum S = 1-spin 
operator placed at each square-lattice point i (Fig. [1]). 
The summations, Z]<(ij))i ^[ijk]^ 

possible rectangular (nearest-neighbor) edges (ij), skew- 
diagonal pairs {{ij)), and skew-diagonal adjacent three 
sites [ijk], respectively. Correspondingly, the coupling con¬ 
stants, J, J', and J", denote the nearest-neighbor-, skew- 
diagonal- and skew-diagonal-adjacent-three-spin-interaction 
parameters. Hereafter, the coupling constant J' is consid¬ 
ered as the unit of energy (J' = 1). The underlying physics 
behind each interaction parameter is as follows. For suffi¬ 
ciently large nearest-neighbor interaction J(> 0), the sys¬ 
tem reduces to a two-dimensional model, and the nematic 
phase emerges [55] ; here, the quadratic component of the 
Heisenberg interaction is set to j = 0.5 throughout this 
study. On the contrary, the coupling constants J'{= 1) 
m and J"{> 0) [551133] strengthen the spatial anisotropy, 
promoting the formation of singlet dimers along the skew- 
diagonal bonds; a schematic phase diagram is presented in 
Fig. 12 In this paper, we incorporate the latter interaction 
term and adjust this (redundant) interaction parameter 
J" as well as the three-spin-interaction component js so 
as to optimize the finite-size-scaling behavior. 
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It has to be mentioned that in the pioneering study 
fM\ . the bipartite-lattice version (without the diagonal 
interaction) of Eq. ([IJ was investigated by means of the 
quantum Monte Carlo method; for the bipartite lattice, 
the spatial anisotropy J' inevitably violates the symme¬ 
try between the horizontal and vertical directions, and the 
asymmetry might alter the nature of the transition |15] . 
Our non-bipartite-lattice version o retains the symme¬ 
try (between the horizontal and vertical directions), as 
would be apparent from Fig. [T] In order to cope with the 
non-bipartite-type lattice, we employ the exact diagonal- 
ization method with the aid of Novotny’s method (screw¬ 
boundary condition) [3T] , which enables us to treat a vari¬ 
ety of system sizes N = 10,12,... in a systematic manner; 
note that the number of spins constituting a rectangular 
cluster is restricted within iV = 4,9,.... 

The rest of this paper is organized as follows. In Sec.lU 
we present the simulation results. The technical details as 
to the screw-boundary condition are presented as well. In 
Sec. [21 we address the summary and discussions, providing 
a brief overview on the past studies of the correlation- 
length critical exponent v. 

2 Numerical results 

In this section, we present the simulation results. To be¬ 
gin with, we explain the simulation scheme to implement 
the screw-boundary condition, namely, Novotny’s method 
m, briefly. Owing to this scheme, we are able to treat 
an arbitrary number of spins, N = 10,12,..., constituting 
a two-dimensional cluster. The linear dimension L of the 


cluster is given by L = y/N, because the N spins form a 
rectangular cluster. 

The Hamiltonian (HD has been investigated extensively 
in some limiting cases. In order to elucidate the phase di¬ 
agram, Fig. 121 we recollect a number of related studies 
we also address a brief account of the pa¬ 
rameter range surveyed in our preliminary study. In the 
limit J —>■ oo, the model reduces to the two-dimensional 
Heisenberg model with the biquadratic interaction. Ac¬ 
cording to Ref. [22], around j = 0.5, the nematic phase 
is realized. With J = 0 and J" — 0, the one-dimensional 
biquadratic-interaction Heisenberg model is realized, and 
the VBS phase emerges m- Similarly, with J = 0 and 
J' = 0, the VBS phase is realized for sufficiently large 
j3 > 0.111 m- Hence, the interaction parameter J in¬ 
terpolates smoothly these limiting cases, and the phase 
diagram, Fig. [21 follows. In the preliminary stage of the 
research, we dwelt on the subspace J" = 0, which was 
studied in Ref. m- Turning on the interaction J" grad¬ 
ually, we arrive at the optimal regime J" « 0.08, as indi¬ 
cated in Fig 121 At least for a moderate range of J" < 0.2, 
a fundamental feature of the phase diagram, Fig. [21 is kept 
maintained; namely, no signature such as an appearance 
of a certain intermediate phase could be detected. 

2.1 Simulation method: Screw-boundary condition 

In this section, we explain Novotny’s method m to imple¬ 
ment the screw-boundary condition. The screw-boundary 
condition enables us to treat a variety of system sizes 
N = 10,12,...; note that naively, the system size is re- 
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stricted within quadratic numbers N = 4,9,... for a rect¬ 
angular cluster. 

In this paper, we follow the simulation scheme reported 
in Ref. [5S] , where the J" term of the Hamiltonian m was 
not taken into account. The missing term is incorporated 
by the addition of the following term to Eq. (5) of Ref. 

ED; 

N 

j" • S,+1 + J3[(S. • S,+i)(S,+i • S,+2) + h.c]}. ( 2 ) 

i=l 

(The index i runs over a one-dimensional alignment i = 
1 , 2 ,..., Ai in a way intrinsic to the screw-boundary con¬ 
dition.) Thereby, we diagonalize the Hamiltonian matrix 
given by Eq. (5) of Ref. [55] with the term ([5|), employ¬ 
ing the Lanczos algorithm for a finite-size cluster with 
N < 20 spins. Rather technically, the diagonalization was 
performed within the zero-momentum subspace, at which 
the magnetic- (triplet-) excitation gap AE opens. 

2.2 Finite-size scaling of AE: Critical point 

In this section, based on the simulation method explained 
in Sec. o we evaluate the excitation gap AE. Thereby, 
we estimate the location of the critical point via the scaling 
analysis oi AE. 

In Fig. 131 we plot the scaled energy gap LAE for 
various J{/J') and N = 10,12,..., 20 {L = '/N); here, 
the interaction parameters are set to J"{/ J') = 0.08 and 
ja = 0.45. According to the finite-size-scaling theory, the 
intersection point of LAE indicates a location of the crit¬ 
ical point Jc « 0.3, because the scaled energy gap should 
be scale-invariant (dimensionless) at the critical point. 


In Fig. 01 we plot the approximate critical point Jc(Li, L 2 ) 
for [ 2 /(Li -h L 2 )f with 10 < A^i < A ^2 < 20 (Li ,2 = 
y^A^i. 2 ); the interaction parameters are the same as those 
of Fig. 131 Here, the approximate critical point Jc(Ti,T 2 ) 
is defined by the formula 

LiAE{Li)\j^j^(^l^^l^) = L2AE{L2)\j=j^(Li,L2)^ ( 3 ) 

for a pair of system sizes {Li., L 2 ). The least-squares fit to 
the data in Fig. |4] yields Jc = 0.2998(45) in the thermo¬ 
dynamic limit L —>■ 00 . This extrapolated value does not 
affect the subsequent analysis (Sec. 12.31) . and we do not 
go into the discussion of the extrapolation error; actually, 
the approximate critical point Jc{Li,L 2 ) (rather than the 
extrapolated Jc) is fed into the formula, Eq. O- 

Last, we address a number of remarks. First, in Fig. 01 
the finite-size behavior seems to be oscillatory; actually, 
the data exhibit successive bumps for quadratic values of 
N{= L?) Ri 9,16,.... Such an oscillatory deviation is an 
artifact m of the screw-boundary condition. Second, the 
set of the coupling constants J" = 0.08 and ja = 0.45 
were determined so as to optimize the finite-size behavior 
of Fig. 0] (as well as Fig.[3]mentioned afterward). Actually, 
as for J" 0.08 and J 3 yf 0.45, the finite-size behavior of 
Jc{Li,L 2 ) suffers from steep finite-size drift and enhanced 
bumps. Last, in the scaling analysis, Fig. |31 we postulated 
the dynamical critical exponent z = 1. Here, we followed 
the conclusion z = 1 obtained in the pioneering studies 01 
0 ]. 
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2.3 Correlation-length critical exponent v 

In this section, we estimate the correlation-length criti¬ 
cal exponent v. Based on the approximate critical point 
Jc{Li,L 2 ) ([3]), we are able to calculate the approximate 
correlation-length critical exponent 

^n{Li/L2) 

" \n{dj[L,AEiL,)]/dj[L2AE{L2)]}\j=j^^L,,L,)^ 

(4) 

for a pair of system sizes {Li,L 2 ). In Fig.jSJ as the symbol 
-h, we plot i'{Li,L 2 ) for [2/(Li -|- ^2)]^ with 10 < < 

N 2 < 20 (Li ,2 = •\/A^i, 2); here, the interaction param¬ 
eters are the same as those of Fig. [31 The data exhibit 
an oscillatory deviation (bump at iV Ri 16), which is an 
artifact of the screw-boundary condition [31]. The least- 
squares fit to the data in Fig. 0 yields u = 0.889(10) in 
the thermodynamic limit. 

Similar analyses were carried out independently for 
various values of J" with ja = 0.45 fixed. As a conse¬ 
quence, the approximate critical point is plotted in Fig. 0 
for (x) J" = 0, {*) 0.04, and (□) 0.12. The least-squares 
fit to these data yields u = 0.8642(87), 0.8823(93), and 
0.883(12), respectively, in the thermodynamic limit. No¬ 
tably enough, the interaction parameter J" governs the 
convergence of v{Li, L 2 ) to the thermodynamic limit. Par¬ 
ticularly, in the optimal regime 0.04 < J" < 0.12, the 
extrapolation to the thermodynamic limit can be taken 
reliability, allowing us to appreciate the critical exponent 
V K, 0.88 rather systematically. For exceedingly large J" > 
0 .12, there emerge a notable finite-size drift (to the counter 
direction) as well as enhanced oscillatory deviations (steep 


bumps around N k, 16) as to v{Li, L 2 ) and even Jc(Li, L2), 
preventing us from extrapolating v unambiguously. Sur¬ 
veying various J" as well as js, we observed that the un¬ 
certainty of the extrapolation is bounded by Av = 0.03. 
As a result, we estimate the correlation-length critical ex¬ 
ponent as 

i/ = 0.88(3). (5) 

The result is consistent with the estimate ly = 0.92(10) 
reported in Ref. |26] , where the three-spin interaction was 
not taken into account. 

In Fig. [S] we assumed implicitly that the dominant 
scaling correction should obey the power law 1/L^. How¬ 
ever, as mentioned in Introduction, there might be a loga¬ 
rithmic correction [221123] , which is not negligible for large 
system sizes. We cannot exclude a possibility that such a 
correction gives rise to an ambiguity as to the estimate 
®. Here, aiming to provide a crosscheck, we carry out an 
alternative analysis of criticality as presented in the next 
section. 

2.4 Scaling plot for AE 

In this section, we present the scaling plot for AE as a 
cross-check of the analyses in Sec. I2.2l and l2.3l 
In Fig. 13] we present the scaling plot, (J — 

LAE for the same interaction parameters as those of Fig. 
|3| Here, the scaling parameters are set to Jc = 0.2998 (Sec. 
12 .21) and i/ = 0.88 (Sec. 12.31) . The scaled data appear to fall 
into a scaling curve satisfactorily for an appreciable range 
of J, confirming the validity of the analyses in Sec. l2.2l and 
nn As mentioned in Introduction, corrections-to-scaling 
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behavior [2211^ may depend on each lattice realization. In 
this respect, we stress that the present S = 1-spin model 
is less affected by severe scaling corrections. 

This is a good position to address a remark on the va¬ 
lidity of the deconfined-criticality scenario. As mentioned 
in Introduction, it is still unclear whether the (Neel-VBS) 
phase transition is continuous or not. The scaling analy¬ 
sis in Fig. E] suggests that the (spatial-anisotropy-driven 
nematic-VBS) phase transition would be continuous. A 
notable point is that the estimate z/ = 0.88(3) JS]) takes 
a rather enhanced value; actually, it appears to be larger 
than, for instance, that of the Heisenberg universality class, 
z/ = 0.7112(5) [35]. Because of the following reason, such 
a feature might exclude a possibility of the discontinuous 
phase transition. As a matter of fact, the discontinuous 
phase transition does exhibit a pseudo-critical behavior 
(for finite system sizes) with an enhanced effective specific- 
heat critical exponent a = du {d: spatial dimension) |33j . 
Resorting to the hyper-scaling relation a = 2 — 3z^, one 
arrives at a suppressed v = 0.4, which is inaccordant with 
ours. Such a tendency is reasonable, because the discon¬ 
tinuous transition is accompanied with a latent-heat re¬ 
lease, enhancing the specific-heat critical exponent a to 
a considerable extent. A comparison with other existing 
estimates for u is presented in the next section. 

3 Summary and discussions 

The deconfined criticality between the nematic and VBS 
phases for the two-dimensional S = 1-spin model with 
both three-spin and biquadratic interactions o was in¬ 


vestigated with the numerical diagonalization method for 
a finite-size cluster with N < 20 spins; so far, the case 
without the three-spin interaction has been investigated 
extensively [21126]. The extended interaction-parameter 
space enables us to search for a regime of suppressed finite- 
size corrections. Actually, the interaction parameter J" 
governs the convergence to the thermodynamic limit, as 
shown in Fig. [1 In fact, for 0.04 < J" < 0.12, the extrap¬ 
olation can be taken reliably. As a result, we estimate the 
correlation-length critical exponent as v = 0.88(3); the 
result agrees with the preceeding estimate v = 0.92(10) 

m- 

As a comparison, we recollect related studies of the 
deconfined criticality, placing an emphasis on the criti¬ 
cal exponent v. First, as mentioned in Introduction, the 
J-Q model has been investigated extensively. Pioneering 
studies reported v = 0.78(3) [1 and 0.68(4) [1. Possi¬ 
bly, there emerges the logarithmic corrections to scaling 
in the vicinity of the deconfined criticality mm , pre¬ 
venting us from identifying the character of the phase 
transition definitely. A recent simulation result for the 
honeycomb-lattice model indicated a rather suppressed 
value of z/ = 0.54(5) [T3]. Second, a unique approach to the 
deconfined criticality was made by the extention of the in¬ 
ternal symmetry to SU(V) with the continuously variable 
N [31]. There was reported a clear evidence of the phase 
transition of a continuous character with the enhanced 
critical exponent v = 0.75-1. Third, the “fermionic” crit¬ 
icality [ISlfTO] indicates v = 0.80(3) [20] and ~ 0.88 [2T] . 
Last, we overview simulation results for the classical coun- 
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terparts. For the dimer model, the critical exponent was J 

J 

estimated as z/ = 0.73(5) [3S], 0.5 [3B], 0.5-0.53 [37], and 

_J 

/i\y' 


Y- _^ 

/ / 

/ / 

A y 

/ 

0.5-0.6 [33]. The hedgehog-suppressed 0(3) model indi- 


/ ^ 

[ijy^ 


cates an enhanced exponent v = 1.0(2) [39]. These con- 






elusions have not been settled yet. Nonetheless, our result 


indicates a tendency toward an enhancement of v such as Fig- 1- The square-lattice model with the skew-diagonal in- 

the fermionic criticality [2011^ . teraction is considered; see the Hamiltonian JT]). The nearest- 

The present simulation result supports the deconfined- neighbor interaction J stabilizes the nematic phase m 

...... • j. 1 j. r J .1 / i- 1 • j. j • \whereas the skew-diagonal biquadratic and three-spin inter- 

criticahty scenario at least tor the (spatial-amsotropy-driven) & e 


nematic-VBS phase transition; see the discussion in Sec. 
I2.4l as well. In Ref. [40], a novel two-dimensional quantum- 
spin model with the three-spin-interaction was introduced. 
The model may be a good clue to the realization of the 
deconfined criticality without the spatial anisotropy. How¬ 
ever, its rich phase diagram has not been fully clarified. A 
progress toward this direction is remained for the future 
study. 
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